Chaotic Quantization: 

Maybe the Lord plays dice, after all? 
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Abstract. We argue that the quantized non-Abehan gauge theory can be obtained as 
the infrared hmit of the corresponding classical gauge theory in a higher dimension. 
^ I We show how the transformation from classical to quantum field theory emerges, and 

calculate Planck's constant from quantities defined in the underlying classical gauge 
theory. 



1 Introduction 



> 

The question, how gravitation and quantum mechanics can be merged into a 
CO , consistent unified theory of all fundamental interactions, is still open. Logically, 

either general relativity (GR), or quantum mechanics (QM), or possibly both, 
will have to be replaced by a different theory at a more fundamental level. The 
' almost universally accepted notion is that it is GR which needs to be replaced, 

. while QM presumably provides a truly fundamental description of nature. Super- 

' string theory, describing our four-dimensional space-time as the low-energy limit 

of a ten- or eleven-dimensional theory, is widely accepted as the most promis- 
ing approach, but neither the precise form nor the full content of this theory is 
, entirely understood at the present time. 

' What about the other option, considering QM as the low-energy limit of a 

more fundamental theory? This question has been raised by 't Hooft, who con- 
jectured that quantum mechanics can logically arise as the low-energy limit of 
a microscopically deterministic, but dissipative theory [1,2]. Explicit, but highly 
' simplified examples for such a mechanism have been constructed [3,4] . In a recent 

publication [5] with Biro and Matinyan, we showed how (Euclidean) quantum 
field theory can emerge in the infrared limit of a higher-dimensional, nonlinear 
classical field theory (Yang-Mills theory). We called this phenomenon chaotic 
quantization to distinguish it from the formal technique named stochastic quan- 
tization [6], not realizing that this term was already introduced by C. Beck 
several years earlier [7] for essentially the same mechanism. What is special 
about Yang-Mills fields, however, is that they "quantize themselves" , as we shall 
discuss below. In Sect. 2, we introduce the concept of chaotic quantization of a 
system with one degree of freedom, which we extend to field theory in Sect. 3. 
In Sect. 4 we review the chaotic properties of classical Yang-Mills theory, be- 
fore we analyze their chaotic self-quantization in Sect. 5. In the final section, we 
enumerate and discuss several open problems. 
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2 Chaotic quantization 

A classical physical system encodes much more information than the analogous 
quantized system. Consider, for instance a point particle in one dimension. The 
classical system is defined by the pair of coordinates {x,p), implying that every 
point in the continuous phase space represents a different state. For the quantum 
system, on the other hand, the uncertainty relation limits the localization of the 
state in phase space to the finite element AxAp ^ h. This observation suggests 
that it may be useful to consider classical systems, whose internal dynamics 
results in a self-afHicted loss of information. 

Deterministically chaotic systems satisfy this condition. For such a system, 
the rate of information loss is encoded in the Lyapunov exponent A, defined as 
the loarithmic rate of divergence between neighboring trajectories: 

\x,{t) - X2{t)\ e^' (A>0) (1) 

or equivalently in the eigenvalue 7 of the Perron-Frobenius operator, defined as 
the logarithmic rate of convergence of the phase space density to its stationary 
limit: 

pix, t) ^ pM^{x) + p'(x)c-T* (7 > 0) . (2) 

Before we pursue this idea further, let us recall the method of stochastic 
quantization [6]. Consider a quantum field (l){x) in Euclidean space with action 
5* [(/)]. The domain of 4> is formally extended into a fifth dimension denoted by r. 
If the field 0(a;, r) obeys the Langevin-type equation 

d SS 

— (t){x,T)=^ — {x,T)+^{x,T), (3) 

where ^(x, r) represents local white noise defined by the moments 

{aX,T))=0, {aX,T)ax',T'))=2S{x-x')S{T-T'), (4) 

then the long-time average of any physical observable converges to the quantum 
mechanical vacuum expectation value: 

hm 1 r dTO[ct>{x, r)] = (O[0(x)])qm • (5) 

T->oo 1 Jq 

Beck's [8] suggestion was to replace the artificial white noise ^ with the "noise" 
generated by a deterministic, but chaotic (more precisely: (^-mixing [9]) process. 

Following Beck [10,11], let us start by considering a dynamical system with 
two variables x,y, which evolves in discrete time steps of length r. We denote 
the state of the system at t„ = nr as (a;„, y„). We are interested in the dynamics 
in the "physical" variable y, if the motion in x is chaotic on short time scales. 
We define the evolution of the system as follows: 



{Xn+l,yn+l) = f{Xn,yn) = {T{Xn),Xyn+T ' Xn) ■ 



(6) 
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Here the map T is assumed to be (^-mixing [9] . Equation (6) can be considered 
as the stroboscobic map of the differential equation 

oo 

y = -iy + T^''^^Xn-i5{t-nT), (7) 

n=l 

with Xn+1 = T{xn) and A = cxp(— 7r). Obviously, the variables {xn} take the 
role of the noise in this equation. 

The Langevin equation (7) is equivalent to the Perron-Frobenius equation 
for the evolution of the phase space density of an ensemble of systems: 

Pn+,{x,y)= xi9^i= 2. xmd^l 

{x,v)ef-Hx',v') ' ' ' xeT-i(x') ' ' ' 

Expanding (8) into powers of (7t)^/^ = r^/^, taking the limit r ^ 0, and 
interpolating pn{x,y) to a function p{x,y,t) which depends continuously upon 
time, 

p{x, y, t) = (p{x, y, t) + f^/'^a{x, y, t) + fb{x, y, t) + f^/'^c{x, y,t) + -- - , (9) 
yields a set of coupled equations for the coefficient functions: 

xeT-i(x') ' ' ' 
a{x',y,t) = pTs^ ("(^'^'*)~^^^(^'^'*)) 

etc. (10) 

Being interested in the dynamics of the physical variable y only, we define 
projected functions 



Po{y,t) 



J dx(p{x,y,t); a{y,t) = J dxa{x,y,t) ; etc. (11) 



For complete maps T, it is possible to show that 

xeT-i(x') ' ' ' 

for all y and t implies 

dxf{x,y,t)= I dxg{x,y,t). (13) 



The first equation in the expansion of the Perron-Frobenius equation in powers 
of f^/^ then becomes a tautology, while the second one takes the form 



dyj 



dxxip{x,y,t) = 0. (14) 
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The desired dynamical equation for the phase space density po{y,t) is obtained 
as the third equation (at order f): 

— J dxxa{x,y,t) = —{ypoiy,t)) + --^ J dxx^ip{x,y,t) - . (15) 

If h{x) is an invariant of the map T, (f{x,y,t) = h{x)po{y,t) is a solution of 
the first of the set of equations (10) for any function Po{y,t), and (14) requires 
(a;) = J dxxh{x) = 0. Finally, (15) turns into 

— Po(2/,t) = ^ (yPo(2/,0) + ^-^Po(2/,t) - dxxa{x,y,t). (16) 

The left-hand side and the first two terms on the right-hand side have the form of 
a Fokker-Planck (FP) equation; the last term on the right-hand side represents 
a source term. For symmetric maps, T(x) = T{—x), it is easy to sec that the 
source term vanishes, and the projected phase space density poiy,t) satisfies a 
homogeneous FP equation. The first correction in f^/^, a{y, t), also satisfies a FP 
equation, but in this case, the source term docs not vanish. In other words, the 
deviations in the equation for the exact projected phase space density p{y, t) = 
J dx p{x,y,t) from a FP equation are proportional to f^/^. 

An explicit solution can be found for, c. g., the Ulam map T{x) = 1 — 2x^ 
with X e [—1,1]. The stationary solution of the FP equation is: 



Po{y,t)=Po{y)=i^^ j e-^y ; (17) 
and the full solution written as a power series in f is: 



1/2 



p{y,t) =po{y) 



l + fV2(2j/-^y3 + 0(f; 



(18) 



To obtain the Euclidean Schrodinger equation, instead of the FP equation, we 
need to introduce a potential V{y) and rescale the auxiliary variable x according 
to [7] 

T{x) T{x) e^'^^y)/^ 

/ ^_ ^ 1/2 

y„+i = \yn + T^I'^Xn > ?/„+! = Ay„ -|- ( 2 ) Xn , (19) 

where = {x^) = J dxx^h{x). Identifing po with the Euclidean wavefunction 
tp, we find that it satisfies the imaginary-time Schrodinger equation: 



■|^^ + ^(2/) ^{y,t)^SEi^{y,t)=0. (20) 



The correction linear in r^/^ satisfies the same equation, but with an additional 
source term: 

5Ea(,,0=(^)'''^^(.,i). (21) 
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The complete wavefunction is = '0 + f^/^a + which approaches the 
Schrodinger wavefunction in the limit f = 7r ^ 0. 

The important insight to take away from this derivation is that an appro- 
priate chaotic process can serve as the source of the random noise required for 
the stochastic quantization of a dynamical system if the time scale on which 
the chaotic process randomizes is sufficiently short, so that the corrections are 
negligible. 

3 Extension to field theories 

Can this mechanism of quantizing classical systems be generalized to fields ^{x, t) 
with X being a point in three-dimensional space? How this can be done is most 
easily seen, when the field theory is defined on a lattice, rather than a spatial 
continuum. Then all one needs to do is introduce a map T together with some 
internal space {^*} at each lattice point Xi. Beck has proposed to define the 
evolution law including a nearest neighbor coupling [7] 

d 

C+i = (1 - g)T (C) + ^ 5^ (e^" + C"^") , (22) 
which has the continuum limit 

&=T{e:^)+(fH^^\ (23) 

with appropriate coupling constants g, g' ■ We will not follow this route further 
here and refer to Beck's recent monograph [12]. 

A more natural approach consists in identifying the local internal map space 
with a compact Lie group Qx- The simplest realization of this idea is the SU(2) 
gauge theory, i. e., the Yang-Mills field theory. In this case, the internal degrees 
of freedom (color) of the gauge field provide the local space for the chaotic map, 
and the nonlinear dynamics of the gauge field uniquely defines the map. As we 
shall see below, there is no need to introduce new degrees of freedom beyond 
those provided by the gauge field (in one additional dimension) itself [5]. Before 
exploring this idea in more detail, however, we need to review what is known 
about the chaotic dynamics of Yang-Mills fields. 

4 Interlude: Chaotic properties of Yang-Mills fields 

The chaotic nature of classical non-Abelian gauge theories was first recognized 
twenty years ago [13,14]. Over the past decade, extensive numerical solutions of 
spatially varying classical non-Abelian gauge fields on the lattice have revealed 
that the gauge field has positive Lyapunov exponents that grow linearly with the 
energy density of the field configuration and remain well-defined in the limit of 
small lattice spacing a or weak-coupling [15,16]. More recently, numerical stud- 
ies have shown that the (3 -|- l)-dimensional classical non-Abelian lattice gauge 
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theory exhibits global hyperbolicity. This conclusion is based on calculations of 

the complete spectrum of Lyapunov exponents [17] and on the long-time statis- 
tical properties of local fluctualtions of the Kolmogorov-Sinai (KS) entropy in 
the classical SU(2) gauge theory [18]. 

It is useful to note some important relationships between ergodic and pe- 
riodic orbits for globally hyperbolic dynamical systems. The ergodic Lyapunov 
exponents Xr,i are obtained by numerical integration of a randomly chosen er- 
godic trajectory, denoted by its origin r. In a Hamiltonian hyperbolic dynamical 
system with d degrees of freedom the sum of its d — 1 positive ergodic Lyapunov 
exponents is obtained as the ergodic mean of the local expansion rate: 

lim hr{t) = lim - / x(x(t')) dt' = V Xr,i = /iks • (24) 

t— >(X) t— »cx) t In 

Here /iks denotes the Kolomogorov-Sinai entropy and 

x(»W) = i.nde.(|2|) (25) 



dt 



expanding 



is the local rate of expansion along the trajectory x{i}j. Due to the equidistribu- 
tion of periodic orbits in phase space it is possible to evaluate the ergodic mean 
in (24) by weighted sums over periodic orbits p. In fact, for hyperbolic systems 
the thermodynamic formalism allows to express certain invariant measures on 
phase space in terms of averages over periodic orbits [19,20]. Labeling periodic 
orbits by p (rather than a starting point), and denoting their periods and pos- 
itive Lyapunov exponents by Tp and Ap.i, respectively, the connection between 
the positive ergodic Lyapunov exponents and those of periodic orbits is given 
by: 

ftKS = ^lim ^ j—^^ ^ ' 

where e > is a small number. The topological pressure P{P) is a useful tool for 
analyzing invariant measures on phase space in terms of periodic orbits. This 
function can be expressed as 



P{(3)= lim Yl exp(-/3^Ap,.TpJ . (27) 

t<Tp<t+e \ i=l / 

The relation (26) then follows from (24) and from the identity — _P'(1) = /iks- 

In order to apply this formalism to the Yang-Mills field, the gauge theory 
needs to be formulated as a Hamiltonian system on a spatial lattice. How to 
do this is well known: The Hamiltonian lattice gauge theory was formulated 
by Kogut and Susskind [21] in order to study the nonpcrturbative properties 
of nonabelian gauge theories, such as quark confinement. Denoting the lattice 
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spacing by a and the nonabelian coupling constant by g, the Hamiltonian for 
the SU(2) gauge theory is given by 

^ = C E ^^x. + i - tr U,,ij) (28) 

x,i x,ij 

Here the Ux,i S SU(2) are called the link variables at point x in the coordinate 
direction i, and the Ex,i E LSU(2) denote the color-electric field strengths com- 
ponents. Ux,ij denotes the plaquette product of four link variables, starting and 
ending at x and circumscribing the elementary square in the i,j directions. The 
chaotic nature of this theory was demonstrated [15] by numerical simulations, 
and Gong [17] obtained the complete Lyapunov spectrum for lattice volumes 
with L = 1,2,3. Bolte et al. [18] extended these calculations to the lattices of 
size L = 4,6. 



0.3 r 




2 4 6 



Fig. 1. Distribution of numerically obtained ergodic Lyapunov exponents for a classical 

SU(2) gauge theory on lattices of size L = 2.4,6. The index i numbers the Lyapunov 
exponents and the abscissa is scaled with L'^ . The energy per plaquette was chosen as 
1.8/g^a. 

For sufficiently long trajectories and fixed energy per lattice site the Lya- 
punov spectrum shows a unique shape, independent of the lattize size, as shown 
in Fig.l. Indeed, for a completely hyperbolic system, physical intuition dictates 
that the KS entropy /iks is an extensive quantity. For this to be true, the sum 
over all positive Lyapunov exponents must scale like the lattice volume and 
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the shape of the distribution of Lyapunov exponents must be independent of L. 
Figure 1 confirms this expectation. 



150 




-150 I ' ' ' ' ' 1 

2x4 6 



Fig. 2. Top: The distribution of the sum of local expansion rates hr{ts) for L=4 and a 
short sampling time ts = 0.1, together with a Gaussian fit. Bottom: The autocorrelation 
function o(t) for this distribution. 

The KS entropy is a global property of the dynamical system. The next step 
of detail of the ergodic nature of the system is provided by the fluctuations in 
the quasi- local average of x{x{t)). These fluctuations are obtained by integrating 
(25) up to a sampling time tg. For sampling times tg much longer than the largest 
correlation time one expects that obscrvablcs sampled along ergodic trajectories 
exhibit Gaussian fluctuations about their ergodic mean. Waddington [22] has 
shown that for Anosov systems (i.e., fully hyperbolic systems on compact phase 
spaces) the probability distribution for /ir(is) is Gaussian with mean hKS- 

and square variance proportional to P"(l): 



Ahr{t) for t^oo 



(30) 
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The variance (30) can be related to the autocorrelation function 

air) = {x{x{r))xixm)-his (31) 
of the local ergodic Lyapunov exponents, through the relation 

t {Ahr{t)f = p a(r) dT ^ P"(l) . (32) 

Figure 2 (top) shows the distribution of sampled valued of hr{ts), (obtained 
on a single, very long trajectory) for a short sampling interval ts = 0.1 for a 
L = 4 lattice. The bottom part of Fig. 2 shows the autocorrelation function 
a(r) obtained for the same trajectory. A numerical fit of the function indicates 
that a(T) decays exponentially for large times. The value predicted for Ahr by 
(30,32) is in excellent agreement with the value obtained by a Gaussian fit to 
the sampled distribution. 



0.4 




0.1 
0.05 



20 40 60 80 100 

Fig. 3. Ahr{ts), scaled with {ts/L^Y^'^ , as a function of te- 

One can also study how the width of the Gaussian scales with the lattice size 
L. To a very good approximation one finds that it is proportional to If 
one includes the dependence on the sampling time, the variance of /i,- scales like 
^Jl^Jts, (see Fig. 3). As the mean value /iks scales like L^, this result confirms 
the Gaussian nature of the fluctuations: Zi/ir//iKS ~ (L^ig)^^^^- 

These numerically obtained results provide strong (although not mathemat- 
ically conclusive) evidence that the SU(2) lattice gauge theory is a strongly 
chaotic (Anosov, (^-mixing) system with properties required for the formalism 
of Sect. 2. 
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5 Yang-Mills fields quantize themselves 

The results discussed in the previous section imply that correlation functions of 
physical observables decay rapidly, and that long-time averages of observables for 

a single initial gauge field configuration coincide with their microcanonical phase- 
space average, up to Gaussian fluctuations which vanish in the long observation 

— 1/2 

time limit astg ■ Since the relative fluctuations of extensive quantities scale as 
the microcanonical (fixcd-cncrgy) average can be safely replaced by the 
canonical average when the spatial volume probed by the observable becomes 
large. In the following we discuss the hierarchy of time and length scales on 
which this transformation occurs [5]. 

According to the cited results, the classical non-Abelian gauge field self- 
thermalizes on a finite time scale Teq given by the ratio of the equilibrium entropy 
and the KS-entropy, which determines the growth rate of the course-grained 
entropy: 

Teq = "Seq/ZlKS • (33) 

At weak coupling, the KS-entropy for the (3-|-l)-dimensional SU(2) gauge theory 
scales as 

hKS - g^E ^ g^T{L/af , (34) 

where E is the total energy of the field configuration and T is the related tem- 
perature related to E by (for SU(A^c) 

e^^=2{N^-l)^+0{g'). (35) 

The equilibrium entropy of the lattice is independent of the energy and propor- 
tional to the number of degrees of freedom of the lattice: ^eq ~ {L/af. The time 
scale for self-equilibration is thus given by 

^ (5'T)-i . (36) 



When one is interested only in long-term averages of observables, it is thus suffi- 
cient to consider the thermal classical gauge theory on a three-dimensional spa- 
tial lattice. Furthermore, on time scales t ^ Teq, the Yang-Mills field generates 
a random Gaussian process, which is required for chaotic quantization. 

The dynamic properties of thermal non-Abelian gauge fields at such long 
distances have been studied in much detail [23-25]. While these studies have 
been made exclusively for the thermal quantum field theory, their results are 
readily transcribed to the thermal classical gauge field with a lattice cutoff. 
The real-time dynamics of the gauge field at long distances and times can be 
described, at leading order in the coupling constant g, by a Langevin equation 

dA 

a,— = -DxB + ^, (37) 
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where D is the gauge covariant spatial derivative, B = DxA\s the magnetic 
field strength, and ^ denotes Gaussian distributed (white) noise with zero mean 
and variance 

{ii{x, t)ij{x', t')) = 2acTSijS\x - x')S{t - t') . (38) 

Here CTc denotes the color conductivity [26] of the thermal gauge field which is 
determined by the ratio oj^/j of the plasma frequency u and the damping rate 
7 of a thermal gauge field excitation. In the classical field theory with a lattice 
cutoff, the color conductivity scales as 

CTc {aln[djna.g/del])~^ . (39) 

This relation implies that the color conductivity is an ultraviolet sensitive quan- 
tity, which depends on the lattice cutoff a. 

We now consider observers measuring physical quantities on long time and 
distance scales {t, a, {g'^T)~^). The random process defined by the Langevin 

equation (37) generates three-dimensional field configurations with a probability 
distribution P[A] determined by the Fokker-Planck equation 

where W[A] denotes the magnetic energy functional 

W[A] = J d^x ^Bix)^ . (41) 

Any non-static excitations of the magnetic sector of the gauge field, i.e. mag- 
netic fields B{k) not satisfying k x B = 0, die away rapidly on a time scale 
of order (Xc/fc^, where k denotes the wave vector of the field excitation. Long- 
term averages are determined by the static magnetic field sector weighted by the 
stationary solution of the FP equation (40): 

Po[A] =e-^[^l/^. (42) 

The observer measures 

{0[A]) = J VAO[A]e-^^^^/^ . (43) 

The magnetic field 

Bi = leijkF^" ^ ^ei^kP" (44) 

defines a three dimensional field strength tensor f^'^, and W/T can be identified 
with the three-dimensional action measured in units of Planck's constant 
n = aT [5]: 
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The rescaling of the gauge field strength by the fundamental length scale a is 
required for dimensional reasons. The same rescaling also determines the three- 
dimensional coupling constant to be 

The central result of this section is that the highly excited classical (3 + 1)- 
dimensional Yang-Mills theory reduces to a vacuum quantum Yang-Mills theory 

in three Euclidean dimensions for an observer who is only interested in physics 
at long distance and time scales, with vacuum expectation values of the standard 
form: 

{0[A])3= j VAO[A]e-^^^^'^'^ . (47) 

Planck's constant is determined by two microscopic quantities of the "funda- 
mental" theory, a and T: 

h = aT. (48) 

The existence of both, a fundamental length scale and an energy scale, is critical 
to the emergence of a constant with the dimensions of an action. 
Some noteworthy comments: 

• It is important to note that the dimensional reduction is not induced by a 
compactification of the time coordinate, cither in real or imaginary time. We 
have not asumed a thermal ensemble of gauge fields in the original Minkowski 
space theory, and the random solution of the (3 + l)-dimensional classical 
field theory does not satisfy periodic boundary conditions in imaginary time. 
The effective dimensional reduction is not caused by a discreteness of the 
excitations with respect to the time-like dimension, but by the dissipative 
nature of the (3 + l)-dimensional dynamics. Magnetic field configurations 
satisfying D x B = can be thought of as low-dimensional attractors of the 
dissipative motion, and the chaotic dynamical fluctuations of the gauge flcld 
around the attractor can be consistently interpreted as quantum fluctuations 
of a vacuum gauge field in 3-dimensional Euclidean space. 

• The dimensional reduction by chaotic; fluc;tiiations and dissipation docs not 
occur in scalar field theories, because there is no dynamical sector that sur- 
vives long-time averaging. Quasi-thermal fluctuations generate a dynamical 
mass for the scalar field and thus eliminate any arbitraily slow field modes. 
An exception may be the case where the excitation energy of the scalar field 
is just right to put the quasi-thermal field at the critical temperature of a 
second-order phase transition, where arbitraily slow modes exist as fluctu- 
ations of the order parameter. In the case of gauge fields, the transverse 
magnetic sector is protected by the gauge symmetry, and it is this sector 
which survives the time average without any need for fine-tuning of the mi- 
croscopic theory. 

• Generalizing the results of Sect. 2 and 4, we expect corrections to the vacuum 
quantum field theory in three dimensions to be of the order of the relative 
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fluctuations of the KS entropy within a spatial region of size Ax: 

^ ~ (g^TAx)-^ ~ ( ^1—^ ' . (49) 
hKS ^ ' WhAx) ^ ' 

If, e.g., a is of the order of the Planck length and Ax is any physically 

accessible length scale, the corrections to the dynamics of the quantized 
Yang-Mills field will be exceedingly small, indeed. 



6 Open problems 

Our example for the chaotic quantization of a three-dimensional gauge theory 
in Euclidean space raises a number of questions: 

1. Does the principle of chaotic quantization generalize to higher dimensions, 
in particular, to quantization in four dimensions? 

2. Can the method be extended to describe field quantization in Minkowski 
space? 

3. Can gravity be included in this framework? Does the nonlinearity of clas- 
sical general relativity provide a source of random noise at short distances, 
allowing for an effective quantum theory to emerge on long distance and 
time scales? 

The first question is most easily addressed. As long as globally hyperbolic 
classical field theories can be identified in higher dimensions, our proposed mech- 
anism should apply. Although we do not know of any systematic study of di- 
cretized field theories in higher dimensions, a plausibility argument can be made 
that Yang-Mills fields exhibit chaos in (4-fl) dimensions. For this purpose, we 
consider the infrared limit of a spatially constant gauge potential [13,14]. For the 
SU(A^c) gauge field in dimensions in the Aq = gauge, there are 3(A'^^ — 1) 

interacting components of the vector potential and 3(7V^ — 1) canonically conju- 
gate momenta, which depend only on the time coordinate. The remaining gauge 
transformations and Gauss' law allow to eliminate 2(Af^ — 1) degrees of freedom. 
Next, rotational invariance in D dimensions permits to reduce the number of 
dynamical degrees of freedom by twice the number of generators of the group 
S0(£'), i.e. by D{D - 1). This leaves a (£) - l)(2Af2 _ 2 - D)-dimensional phase 
space of the dynamical degrees of freedom and their conjugate momenta. For 
the dynamics to be chaotic, this number must be at least three. For the simplest 
gauge group SU(2), this condition permits infrared chaos in 2 < JD < 5 dimen- 
sions, including the interesting case D = 4. Higher gauge groups are needed to 
extend the chaotic quantization scheme to gauge fields in D > 5 dimensions. Of 
course, this reasoning does not establish full chaoticity of the Yang-Mills field 
in these higher dimensions, it only indicates the possibility. Numerical studies 
will be required to establish the presence of strong chaos in these classical field 
theories. 

The second question is more difficult. A formal answer could be that the 
Minkowski-space quantum field theory can (and even must) be obtained by 
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analytic continuation from the Euclidean field theory. Any observable in the 
Minkowski space theory that can be expressed as a vacuum expectation value of 
field operators can be obtained in this manner. If this argument appears some- 
what unphysical, one might consider a completely different approach, beginning 
with a chaotic classical field theory defined in (3+2) dimensions. Field theories 
defined in spaces with two time-like dimensions were first proposed by Dirac in 
the context of conformal field theory [27] and have recently been considered as 
generalizations of superstring theory [28]. In this case, one time dimension is 
effectively eliminated by gauge fixing. 

In the absence of similar explicit constraints, field theories with two time- 
like dimensions, even if the second time direction is compact, exhibit unphysical 
properties, such as a lack of causality and unitarity [29]. E. g., the Coulomb po- 
tential of a point charge in the presence of a second, curled up timelike dimension 
with period L is complex: 



The lack of causality is closely related to the problem of the existence of time- 
like closed loops, which confuse the distinction between past and future. These 
problems are avoided, if the second timelike dimension is "thermal" , i. e. if it is 
compact in the imaginary time direction. The factor i then disappears from the 
exponent of (50), and the corrections to the Coulomb law are real and exponen- 
tially suppressed at large distances. In the context of the mechanism discussed in 
Sect. 5, the physical time dimension may be defined as the coordinate orthogonal 
to the total 5-momentum vector P'^ of the initial field configuration. Whether 
this reasoning applies to the case, where the "thermal" field theory is really an 
ergodic one, remains to be confirmed, but it is quite plausible. 

In the presence of two time-like dimensions, the "energy" becomes a two- 
component vector E, which is a part of the (_D+2)-dimensional energy- momentum 
vector. If we select an initial field configuration with energy En, where n is a 
two-dimensional unit vector, this choice defines a preferred time-like direction 
tn, in which the field thermalizes. Conservation of the energy-momentum vector 
ensures that the total energy component orthogonal to n always remains zero. 
The choice of an initial field configuration corresponds to a spontaneous breaking 
of the global S0(D,2) symmetry down to a global S0(D,1) symmetry. Whether 
this process leads to an effective quantum field theory in the (Z)-|- l)-dimensional 
Minkowski space, remains to be investigated. 

Finally, what about gravity? One reason, why this question is difficult to 
answer, is that little is known about the properties of general relativity as a 
dynamical system. Due to its different gauge group structure, gravity has more 
"capacity to resist" chaos than the Yang-Mills fields. Local invariance against 
coordinate transformations is incompatible with the concept of Hamiltonian dy- 
namics, which has a preferred time direction. The Hamiltonian version of general 
relativity [30] used in cosmology (mixmaster universe) reflects this exceptional 
situation in GR. Even the most basic definitions of deterministic chaos are not 




(50) 
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directly applicable to general relativity and require appropriate generalizations. 

Many special, chaotic solutions of Einstein's equations have been found [31,32]. 
Most important among these is the eternally oscillating chaotic behavior dis- 
covered in [33] for the generic solution of the vacuum Einstein equations in the 
vicinity of a spacelike singularity, which has the character of deterministic chaos. 
However, it is not known whether the generic solution exhibits chaos, as in the 
case of the Yang-Mills theory. It is not even clear what the proper framework 
for a systematic numerical study of this question would be. For example, the 
Lyapunov exponents, which were so effectively used in Yang-Mills theory, are 
coordinate dependent in GR due to the general covariance of Einstein's equa- 
tions against coordinate transformations. 

What is clear, is that general relativity shares many of the properties, which 
allow nonabelian gauge theories to chaotically quantize themselves: Einstein's 
equations are strongly nonlinear and have a large set of gauge invariances which 
could guarantee the survival of a dynamical sector at long distances in the pres- 
ence of quasi-thermal noise. Under such conditions, GR may not even require 
a short-distance cutoff, because the thermal Schwarzschild radius (2GT) defines 
an effective limit to short-distance dynamics which can couple to the dynamics 
at large distance scales. Applying the relation (48) determining h for the Yang- 
Mills field, one might conjecture that the analogous relation for gravity has the 
form h ~ GT^. If the temperature parameter T were of the order of the Planck 
mass, this relation would yield the observed magnitude of h. 

However, it is not clear whether the same mechanism - chaos, or exponential 
growth of sensitivity to initial conditions - which causes information loss in the 
dynamics of Yang-Mills fields, must also operate in general relativity, 't Hooft 
has speculated that microscopic black hole formation may be the mechanism 
that causes the loss of information in the case of gravity [1]. Again, a much 
better understanding of the structure of generic solutions of Einstein's equations 
must be a prerequisite to an exploration of these interesting questions. 
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